Multiphoton coherent manipulation in large-spin qubits 
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Large spin Mn 2+ ions (S = 5/2) diluted in a non-magnetic MgO matrix of high crystalline 
symmetry are used to realize a six level system that can be operated by means of multi-photon 
coherent Rabi oscillations. This spin system has a very small anisotropy which can be tuned in- 
situ to reversibly transform the system between harmonic and non-harmonic level configurations. 
Decoherence effects are strongly suppressed as a result of the quasi-isotropic electron interaction with 
the crystal field and with the 55 Mn nuclear spins. These results suggest new ways of manipulating, 
reading and resetting spin quantum states which can be applied to encode a qubit across several 
quantum levels. 

PACS numbers: 03.67.-a,71.70.Ch,75.10.Dg,76.30.Da 
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Quantum algorithms offer the possibility of solving 
certain problems with a tremendous increase of speed 
compared to their classical counterparts Manipu- 
lating quantum bits, however, is a delicate process, and 
photons are an ideal choice as they interact with quan- 
tum systems in predictable ways. As qubit implemen- 
tations, spin-based systems are nanoscopic, potentially 
operable up to room temperature and less noise sensi- 
tive and benefit of single-spin detection schemes 0, 0, 0- 
When diluted enough to avoid spin-spin interactions, a 
variety of systems show long coherence times, e.g., the 
nitrogen- vacancy (N-V) color centers in diamonds [a, 0], 
N atoms in CfinUj, Ho 3+ and Cr 5+ ions[^,[I3| and molec- 
ular magnets [111. Il2j. In this Letter, we present first ob- 
servations of multi-photon spin coherent manipulation in 
a multi- level system (Mn 2+ ) diluted in MgO, a highly 
symmetric nonmagnetic matrix. 

Such multi-state systems are proposed for quantum 



algorithms in either size-limited [13| or scalable [14[ 
schemes and to study new exoti c q uantum phenomena, 
like the quantum antiresonances [1 51 ] , by using microwave 
(MW) pulses to generate entangled states. The high sym- 
metry of the host ensures that the Mn 2+ spins are in an 
almost isotropic environment, which permits reaching a 
regime where the MW excitation [last term of Hamilto- 
nian |[T])] and the induced Rabi splittings are comparable 
to spin's anisotropy. Mn 2+ ions are diluted in a single 
crystal of MgO with cubic symmetry F m ^ m (lattice con- 
stant 4.216 A). They replace Mg in a concentration esti- 
mated at 1 ppm. The spin Hamiltonian is given by (refs. 
US): 



H = a/6 [Si + S^ + St- S{S + l)^ 2 - l)/5] (1) 
+gii B H .S - AS. I + gfj, B h mw .S cos(27r/t) 

where S x ^ VyZ are the spin projection operators, g — 2.0025 
is the gyromagnetic ratio, a = 55.7 MHz is the anisotropy 
constant, A — 244 MHz is the hyperfine constant of 
55 Mn (/ = 5/2), h mw and / represent the MW ampli- 



tude and frequency respectively, and Hq is the static field 
(Hq _L h rnw ). The g, A and a parameters were deter- 
mined experimentally by continuous wave spectroscopy 
at room temperature and are in excellent agreement with 
literature values [3, [T3|- All Electron Paramagnetic 
Resonance (EPR) measurements are performed using a 
Bruker Elexsys 680 spectrometer at X-band (~9.6 GHz) 
in both continuous-wave (cw) and pulsed modes. 

Spin dynamics under coherent driving is studied using 
the rotating frame approximation [18] (RFA), which con- 
sists of rotating the Hamiltonian's reference frame with 
a frequency / around the direction of the static field Hq. 
Experimental conditions ensures a < A <C gfisHoS/h = 
f = 9.62 GHz, where h is Planck's constant. This implies 
that (i) -ffo's direction can be approximated as a good 
quantization axis for S and H and (ii) coherent MW driv- 
ing is confined between levels of same nuclear spin pro- 
jection mi (see also refs. [H,[l3]). A numerical diagonal- 
ization of the static Hamiltonian [l^ | is performed, and 
the truncation to the desired mi subspace leads to a di- 
agonal Hamiltonian H mi . The later is used to obtain the 
RFA Hamiltonian H rot (see Eq. (30) in ref. [jp]). The 
resulting time-dependent density matrix p(t, T) , where T 
is temperature, gives the ensemble-averaged magnetiza- 
tion (Sz)(t) = Tr (p(t,T)S z ). Also, the energy diagrams 
of the dressed states (spin-states in interaction with the 
microwave field [2lj ] ) are computed by diagonalizing H rot 
as a function of the detuning SH = H — H resi where H res 
is the field of the central resonance S z = —1/2 «-> +1/2 
(e.g., H res = 329.91 mT; see Fig. EE). 

The corresponding diagrams for two h mw values are 
given in Fig. [IJl as a function of SH for i? ||[110] axis. 
At low MW powers, the dressed state diagram shows 
level crossings where resonances in a typical Zeeman di- 
agram appear (hf is subtracted from consecutive levels). 
Higher MW powers induce splittings of size hFn, where 
Fr is the Rabi frequency of the spin's coherent nutation. 
The participating number of photons is indicated near 
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FIG. 1: (a) Dressed state energy diagrams for two h mw val- 
ues, showing level repulsions (or Rabi splittings) when the 
Zeeman splitting matches the integer number of photons la- 
beling each resonance, (b) Level diagram indicating the level 
spin projection (on the left) involved in each multi-photon 
transition, (c) CW absorption derivative at room tempera- 
ture for increasing MW power (from 3 to 23 dBm; curves are 
shifted vertically). At low powers, only the central (1) and 
the laterally detuned single-photon resonances (l',l") are ob- 
served. At high powers, multi-photon resonances are revealed 
(2,2'-unresolved and 3 at the central resonance). 



each resonance and the initial and final S z projections 
are shown in Fig. [TJs. At zero detuning, the one-, three- 
and five-photon resonances are active with the 1-photon 
largely dominating the dynamics [see the size of the Rabi 
splittings in Fig. [1^]. At larger detunings, one observes 
two 1-photon resonances and a 3-photon (3'), a 

two-, a four- and another two-photon (2') resonances of 
relatively low splittings. The existence of quasi-harmonic 
intermediate levels ensures, however, sizeable Rabi split- 
tings even for high numbered resonances. Room tem- 
perature continuous wave (cw) spectroscopy indicates 
the existence of single-, double-, and triple-photon res- 
onances when scanning the detuning of the static field as 
in Fig.QJ. The field scan is performed for the to/ = —3/2 
group of levels. Similar peak structures exist for all other 
To/ projections and are separated by ~8.7 mT (see also 
ref. 



22j). At low MW power, the absorption derivative 
indicates a strong central (1) and two unresolved lateral 
(1' and 1") one-photon resonances. When increasing the 
nominal MW power from 3 to 23 dBm, the signal from 
the single-photon resonances decreases due to a satura- 
tion effect, while double- (2, 2') and triple-photon (3') 
spectroscopy peaks are revealed. The dashed lines indi- 
cate the resonances and are at symmetrical detunings as 
in Fig. [Ik.. 



Coherent Rabi driving 23] is achieved by setting the 
static field at the central resonance and applying one res- 
onant MW pulse. After each pulse, the average magneti- 
zation (S z ) is measured by a spin-echo sequence 12ll24|. 



This ensures an ensemble averaging over the active spin 
packets of the single crystal. (S z ) evolution as a function 
of time is shown in Fig. [2^ for three values of h mw at 
T =30 K. The static field is applied along the unit vec- 
tor (1,0.8, l)/v / 2-64 and from here on data is taken for 
the mi = —3/2 resonances. The Fast Fourier Transform 
(FFT) of the coherent motion indicates peaks for each 
Rabi splitting, as shown in Fig. [2)3 with red dots (one- 
photon) and black squares (three-photon driving). The 
contour plot represents simulated FFT, and an experi- 
mental FFT at the maximum power is shown in Fig. 
The simulations are in excellent agreement with the ex- 
perimental findings. The MW field amplitude is the sin- 
gle fitting parameter giving the h mw values of Fig. [2^. At 
increased powers h mw =1.1 mT and 1.4 mT, the three- 
photon coherent driving is visible as a slow oscillation 
superposed on the fast single-photon one. The 3' Rabi 
splitting is slightly detuned but is visible due to the finite 
linewidth of the resonance. The mixing of several levels 
at high power adversely affects the Rabi decay time, but 
one notes, however, a large amplitude three-photon os- 
cillation. 

The anisotropy term a acts as a small perturbation in 
the spin Hamiltonian but impacts significantly the spin 
dynamics. By means of crystal rotation, one can study 
the role of the anisotropy at fixed static and MW field 
amplitude. Rabi oscillations are recorded each 1° and an- 
alyzed by FFT. The angle dependence of Fig. [3] presents 
experimental FFT peaks in the to/ = —3/2 case (sym- 
bols in black) superimposed on a calculated FFT contour 
plot with h mw as single fit parameter. The temperature 
is T=40 K, and the field direction is along the unit vector 
k — (cos 9, sin 9, 0) [see the inset in Fig.[3k]. As expected 
from the cubic symmetry, the data shows a 90° peri- 
odicity in 9. For \8\ < 20°, the one-photon oscillation 
(squares) dominates the dynamics at a rate ~53.5 
MHz corresponding to h mw =1.3A mT. This identifies the 
dominant transition as being between the S z projections 
— 1/2 — > +1/2 [see resonance (1) in Fig.Q})] in agreement 
with the one-photon Rabi frequency given by (ref. [HI]): 



hF k = -j9VBK 



,^S(S + 1)-S Z (S Z + 1) (2) 



The relatively high value of the MW drive reveals three- 
photon oscillations for 9 closer to 20° and above (the 
black triangles). Of particular interest is the region 
25° < 9 < 65° containing two nodal points at 9 = 31° 
and 59° for which the effect of transverse anisotropy 
is compensated by field's direction [3]. At these com- 
pensation angles anisotropy effects are annihilated and 
the level system becomes harmonic. All levels with 
S z = —5/2 to +5/2 are brought into a coherent superpo- 
sition. Consequently, it behaves as an effective two-level 
system. The measured Rabi oscillation and its FFT at 
the compensation angle [see also the colored markers on 
the ordinate axis of Fig. [3^,] are given in Fig. and c 
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FIG. 2: (a) Rabi oscillations at T=30 K. At low power 
hmm = 0.32 mT, the one-photon oscillations show a decay 
time of ~0.5 /is. At higher powers a slow oscillation is super- 
posed on the one-photon one and is attributed to a coherent 
three-photon oscillation, in agreement with simulations, (b) 
FFT peaks of the experimental Rabi oscillations (red dots 
- single photon, black squares — triple- photon oscillations). 
The full width at half maximum (FWHM) of the FFT peaks 
are the error bars, here similar to the symbol size. The con- 
tour plot shows FFT traces of simulated Rabi oscillations with 
hmm as single fit parameter. The color code goes from (blue) 
to 1 (red) with mixed colors for intermediate values, (c) Ex- 
perimental FFT trace for h mw = 1.4 mT showing two distinct 
peaks. 



respectively, showing a single mode oscillation with fre- 
quency F c ° mv =18.07 MHz. Note that F c ° mv is given 
by the square root's pre-factor of Eq. [2] which depends 
linearly on h mw , in agreement with a two-level picture. 
The experimental ratio F^/ F^ mp = 2.96 is in agreement 
with an expected value of 3 [see Eq.flSJ]. The angular re- 
gion where the anisotropy is compensated is found to 
be narrow, within 1°. For angles slightly shifted away 
from the compensation angle, anisotropy generates a bi- 
modal Rabi oscillation, as the one presented in Fig. [3b 
for 9 = 33°, with a beating period given by modes differ- 
ence. Accordingly, the corresponding FFT indicates two 
peaks located at 17.58 MHz and 22.95 MHz [Fig. [5b]. 
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FIG. 3: (a) FFT peaks, in black, of Rabi oscillations recorded 
at T—40 K, h m w = 1-34 mT, and mi = —3/2, as a function 
of 9 (the inset shows the relative angles between the mag- 
netic fields and the crystal axes). Error bars indicate peaks 
FWHM. The contour plot represents FFT traces of simulated 
Rabi oscillations (color code as in Fig. [2}. The markers on 
the ordinate axis indicate the angles used in (b) and (c). (b) 
At the compensation angle 9 = 31°, the multi-level system 
behaves as a two-level system showing a single-valued Rabi 
frequency. At 9 = 33° beatings between two close Rabi fre- 
quencies are observed. At 9 — 45°, (Sz) undergoes a fast 
nutation superposed on a slower and less defined one. (c) 
FFT of data in (b) [curves are shifted vertically, as in (b)]. 

In the region between the two compensation angles 
31° < 9 < 59°, one observes numerically a complicated 
structure of multi-photon Rabi splittings developing in 
the dressed state energy diagrams. Experimentally, the 
Rabi oscillations do show a multi-mode structure, as the 
one measured for 9 — 45° [Fig. [3b>]. Its FFT indicates 
a narrow peak at 35.7 MHz corresponding to the single 
photon drive and a broad peak centered at ~13 MHz 
comprising three- and five-photon transitions [Fig. [3b] . 

The phenomena described here are visible at room 
temperature (RT) as well, though with less intensitiy. 
As an example, Rabi oscillations at the compensation 
angle 9 = 31° and also at 9 = 45° are given in Fig. 0] for 
the case mi = —1/2. Here, the (S x ) projection is mea- 
sured by the Free Induction Decay (FID) method imme- 
diately after the end of the driving MW pulse [10j. As 
discussed above, the ratio of the Rabi frequencies (here 
21 MHz/7 MHz) fulfills the condition F^/F^ mp = 3. 

Multi-photon oscillations are shorter lived when com- 
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FIG. 4: Rabi oscillations measured at room temperature and 
hmw =0.52 mT, for — 31° (in red) — at the compensation 
angle, when the multi-level systems behaves as a two-level 
system - and for 6 = 45° (#o||[101], in blue). 

pared to the one-photon one due to the larger number 
of levels involved, as visible in Fig. [51 However, the slow 
oscillation attributed to the three-photon resonance re- 
mains coherent up to 0.4 /is. This shows that an im- 
portant part of the noise spectrum is located in the high 
frequency range and affects nutations with Fr > 20 MHz. 
Note that the large hypcrfinc held of the 55 Mn has a very 
limited influence on the decoherence of the electronic 
spin owing to an isotropic hyperfine coupling A. Simi- 
larly, the host's cubic symmetry and Mn 2+ 's zero orbital 
momentum ensure a negligible spin-orbit coupling which 
decreases relaxation processes. Dipolar Mn 2+ — Mn 2+ 
couplings are estimated to be rather weak in this diluted 
system (~ 0.01 MHz) due to an average spin-spin dis- 
tance of ^42 nm. The main source of decoherence resides 
on the superhyperfine coupling between the Mn 2+ spins 
and the 25 Mg nuclear spin bathjIH] (7 = 5/2 with 10% 
natural abundance). 

In this system level coupling and anharmonicity can 
be tuned by field orientation and MW power. For in- 
stance, consecutive MW pulses of different intensities will 
couple different levels via controlled multi-photon pro- 
cesses. Sweeping the field orientation through the com- 
pensation angle where levels are equally spaced could 
provide a means of resetting the qubits or additional gat- 
ing capabilities, depending of the MW power. The use 
of large spins can simplify significantly the implementa- 
tion of quantum algorithms since they provide a natural 
multi-level (or multi-qubit) scheme or can be used to en- 
code qubits across several quantum levels. Coupling can 
be achieved by dipolar interactions at higher concentra- 
tions, by photons when integrated in microcavities, or 
by current oscillations if spin detection is done by elec- 
tronic transport. Large spin qubit readout can be done 
in various ways, for instance as in quantum dotsjH, 0], 
optically [B| or by using surface plasmons 0. Such imple- 
mentations adapted to the case of diluted spins in single 



crystals remain to be explored. 

These first observations of spin states dressed by coher- 
ent photons show the leading role of a small anisotropy on 
the dynamics of a spin qubit. Our results apply to other 
spin qubit implementations with quasi-harmonic energy 
diagram and in quasi-isotropic environments. 
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